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According to the Hugenholtz-Van Hove theorem, nuclear symmetry energy Esym{p) and its slope 
L{p) at an arbitrary density p are determined by the nucleon isovector (symmetry) potential 
Usym (p, k) and its momentum dependence — 4^ . The latter determines uniquely the neutron- 
proton effective k-mass splitting 7n^_p(p, 5) = {m'^ — mp)/m in neutron-rich nucleonic matter of 
isospin asymmetry 5. Using currently available constraints on the Esymipo) and 10 at normal 
density po of nuclear matter from 24 recent analyses of various terrestrial nuclear laboratory exper- 
ts , iments and astrophysical observations, we infer the corresponding neutron-proton effective k-mass 
splitting m*^_p{po,S). While the mid-values of the 7Ti*_p(po, 5) obtained from most of the studies 
^^ ' are remarkably consistent with each other and scatter very closely around an empirical value of 
Cn , ™-n-p(po,'5) = 0.24 • i5, reduced experimental error bars and more complete information about the 
^ ■ uncertainties in extracting the Esymipo) and 10 from data using various models are much needed in 
^L{' order to access more accurately how reliably the inferred m*-^_p{po,5) can be used in investigating 
■^C I many interesting issues in both nuclear physics and astrophysics. 
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I. INTRODUCTION 



The ultimate goal of investigating properties of neutron-rich nucleonic matter through terrestrial nuclear laboratory 
experiments and astrophysical observations is to understand the underlying isospin dependence of strong interaction in 
nuclear medium [1| . The Equation of State (EOS) of neutron-rich nucleonic matter can be written within the parabolic 
approximation in terms of the binding energy per nucleon at density p as E{p, S) = E{p, 6 = 0) -\- i?sym(/o)(5^ + 0{S^) 
where S = (p„ — Pp)/{pp + Pn) is the neutron-proton asymmetry and -Esym(p) is the density-dependent nuclear 
symmetry energy. The latter has important applications in many areas of both nuclear physics, see, e.g., refs. 
f^ ' [2-[a| and astrophysics, see, e.g., refs. |9l4ll|. However, the density dependence of nuclear symmetry energy has 
fT^ , been among the most uncertain properties of neutron-rich nucleonic matter. Predictions using various many-body 
theories and interactions diverge quite broadly especially at abnormal densities. It is thus exciting to see that 
(^ I significant progress has been made recently in constraining the Esym{p) around poi see, e.g., ref. [l2l based on model 
fT^ ■ analyses of experimental and/or observational data. In particular, at least 24 recent studies have extracted the slope 



L(po) = [2ip{dEsyai/dp] and -Esym(po) at po [l3M39l |. It is thus interesting to ask timely what wc can learn about the 
isospin dependence of in- medium nuclear interaction from the extracted constrains on 10 and Esymipo) ■ In this work, 

- ^ , we answer this question at the mean-field level by using a formulism developed earlier based on the Hugenholtz-Van 
1^^ • Hove (HVH) theorem. In particular, we infer both the magnitude of the symmetry potential Usym{po,kp) and the 

2 ' neutron-proton effective k-mass splitting r7T.*_ (po, J) corresponding to each of the 24 constraints on Esym{po) and 

- - - 10 at pq. The consistency of the extracted values for C/sym(poj ^f) and m'^_p{pQ,S) from various constraints is then 

examined. It is found that while the mean values of the Usymipo^kp) and r7T,*_ (po,5) from different studies are 
consistent with each other and most of them scatter closely around Ugymipo, kp) =29 MeV and m* _ (po, (5) = 0.24 ■ S, 
respectively, the individual uncertainties from different analyses are still large. Reducing the experimental error bars 
and better understanding the uncertainties of the model analyses are much needed in order to use reliably the extracted 
mean values of the Usymipoikp) and 77i*_p(po,(5) in solving many important problems in both nuclear physics and 
astrophysics. 
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II. RELATIONSHIP BETWEEN NEUTRON-PROTON EFFECTIVE MASS SPLITTING AND 
SYMMETRY ENERGY BASED ON THE HUGENHOLTZ-VAN HOVE THEOREM 

According to the well-known Lane potential [40|, the neutron/proton (n/p) single-particle potential Un/p{p,k,5) 
can be well approximated by 

[/„/p(p, fc, S) « C/o(p, k) ± Usyn.ip, k)S , (1) 

where the Uo{p,k) and Usyni{p,k) are, respectively, the nucleon isoscalar and isovector (symmetry) potentials for 
nucleons with momentum k in asymmetric nuclear matter of isospin asymmetry 6 at density p. Their momentum 
dependence is normally characterized by the nucleon effective k-mass 

« / r ITT- dUr 111 ,„x 

where r =n, p and for neutrons, protons and nucleons, respectively, and m = {rUn + mp)/2 is the average mass of 
nucleons in free-space. While the nucleon isoscalar potential and its momentum dependence, especially at poi have 
been relatively well determined, our knowledge about the isovector potential Usym{p, k) and its momentum dependence 
— ^^ even at normal density is still very poor. However, from the structure of rare isotopes and mechanism of heavy- 
ion reactions to properties of neutron stars, solutions to many interesting issues depend critically on the nucleon 
isovector potential and its momentum dependence. 

Using the Brueckner theory [41| or the Hugenholtz-Van Hove (HVH) theorem [43], the Esym{p) and L{p) can be 
expressed as [H IH-Iiil 

1 h'^k'^ 1 

^sym(p) = 3-3^ + 2^"^"''-^' ^^-'' '-^'' 

L{P) = l^ + |c/sy„.(p,fcF) + ^\k.k^, (4) 

3 zttiq 2 ok 

where kp = {3n'^p/2y^^ is the nucleon Fermi momentum. We emphasize that these relationships are general and 
independent of the many-body theory and/or interaction used to calculate the Usym {p, k) and m^ . It is well known 
that given both two-body and three-body nuclear interactions, the resulting nucleon potential often depends on the 
many-body theory used. On the other hand, the single-particle mean-field potential is often the one directly tested 
in comparing model calculations with experimental/observational data. For example, it is the input for most shell 
model calculations of nuclear structure and transport model simulations of nuclear reactions. The above expressions 
for Esymip) and L{p) indicate that one can use the density dependence of nuclear symmetry energy extracted from 
experiments/observations to test directly the nuclear isovector potential and its momentum dependence, or vice 
versa, without the hinderance of remaining difficulties and uncertainties in nuclear many-body theories. Here, we 
are interested in learning about the isospin dependence of in-medium nuclear interaction at po from the constrained 
Esymipo) and 10. More explicitly. 



UsymiPO^kp) - 

(Po) 



^^syni 



/ \ i m 
Esymipo) - T,^Ef{Po) 

1 TTl 

L{po) ~ 3Esyni{po) + i^^Epipo) 
3 mS 



(5) 
(6) 



dk 

where Epipo) is the Fermi energy at po. While the Usymipo, kp) is completely determined by Esym{po) and ^, the 
(po) also depends on 10. 



dUs. 



dk 

The nucleon effective mass describes to leading order effects related to the non-locality of the underlying nuclear 
interactions and the Pauli exchange effects in many-fermion systems |46l - |48j . While the nucleon isoscalar effective 
k-mass is well determined to be ttiq/tti = 0.7 ±0.05 at po |47l |. essentially nothing is known about the nucleon isovector 
effective mass J48i | . Knowledge about the neutron-proton effective mass splitting is essential for understanding many 
interesting questions in both nuclear physics and astrophysics pil l49l - [54| . such as, pairing and superfluidity in nuclei, 
properties of rare isotopes, isospin transport in heavy-ion reactions, thermal and transport properties of neutron star 
crust and cooling mechanism of protoncutron stars. Unfortunately, even the sign of the neutron-proton effective mass 
splitting, not to mention its magnitude, has been a longstanding and controversial issue. While some theories predict 
that TO* > m* the opposite has often been shown by studies using different models or interactions, see, e.g., refs. 



TABLE I: Constrained values of Esym {po) and 10 from analyses of terrestrial nuclear experiments and astrophysical observations 

Analysis Esymjpo) 10 Ref. 

Atomic masses and n-skin of Sn isotopes (2011) 

FRDM analysis of atomic masses (2012) 

lAS+LDM analysis of masses (2009) 

Atomic masses and n-skin in an empirical approach (2012) 

IQMD analysis of isospin diffusion at 35 MeV/A (2010) 

IQMD analysis of isospin diffusion at 50 MeV/A (2009) 

IBUU04 analysis of isospin diffusion at 50 MeV/A (2005) 

isoscaling analysis of fragments (2007) 

Liquid drop model analysis of atomic masses (2012) 

Droplet Model+n-skin (2009) 

SHF+n-skin (2010) 

Atomic masses (2010) 

Global nucleon optical potential (2010) 

Pygmy dipole resonances (2007) 

Pygmy dipole resonances (2010) 

Thomas-Fermi model analysis of masses (1996) 

AMD analysis of transverse flow (2010) 

r-mode instability of neutron stars (Vidana, 2012) 

r-mode instability of neutron stars (Wen, 2012) 

Mass-radius of neutron stars- analysis 1 (2010) 

Mass-radius of neutron stars-analysis2 (2012) 

Torsional crust oscillation of neutron stars (Gearheart, 2011) 32.5 ± 7.5 

Torsional crust oscillation of neutron stars (Sotani, 2012) 

Binding energy of neutron stars (2009) 

p, M, Sa-Sa, ISa-lSOl- Thus, a convincing conclusion on this issue will have profound ramifications in both nuclear 
physics and astrophysics. The momentum dependence of the isovector potential is conventionally measured by using 
the neutron-proton effective mass splitting 



30.5 ± 3 


52.5 ± 20 


^ 


32.5 ±0.5 


70 ±15 


M 


32.5 ±1.0 


95 ±19 


[15] 


32.1 


64 ±5 


[16j 


30.1 


52 ±0 


\rn 


32.5 ±2.5 


77.5 ±32.5 


[18, Wj 


30.5 ± 0. 


86 ±25 


[20,^ 


31.6 


65 


[22j 


29.6 ±3. 


46.6 ± 37 


[23j 


31.5 ±3.5 


55 ±25 


[24,^ 


30.5 ±5.5 


41 ±41 


[26] 


31.1 ±1.7 


66 ±13 


[27] 


31.3 ±4.5 


52.7 ±22.5 


[28] 


32 ±1.8 


43 ±15 


[29] 


32 ±1.3 


65 ±16 


[30j 


32.65 


50 


m 


30.5 


65 


[32j 


30. ±5 


> 50 


[33] 


32.5 ±7.5 


<65 


m 


31 ±3 


50 ±10 


[35j 


33 ±1.6 


46 ±10 


[36j 


32.5 ±7.5 


<50 


m 


32.5 ±7.5 


115 ±15 


[38j 


32.5 ±7.5 


< 70 


[391 



m„ 



^{p,,S)^-^ ^-^-2S- 



m dU. 



sym 



According to Eq. |6l the Tn*^_ (po, S) is completely determined by the Esym{po) and 10 via 



, TO , 

(2-r - !)• (7) 



w*_p(Po,<5) ^5- 



1 7Ti 

3£'sym(po) - L{po) - 7:^Ef{po) 



3 ttiq 



Ef{po) ■ (2- 



(8) 



It is clear that whether the m* is equal, larger or smaller than the m* depends on if the condition -L(po) ^ [3£^sym(po) ^ 
Ef{po)\ is satisfied. For example, using the most widely accepted empirical vales of Egymipo) =31 McV, m'Q/m = 



1 m 

3 rriQ 

0.7 MeV and Ef{po) = 36 MeV, to obtain a m*^_p{po, S)>Oa. value of 10< 76 MeV is required. 

III. NEUTRON-PROTON EFFECTIVE MASS SPLITTING FROM CONSTRAINS ON THE DENSITY 
DEPENDENCE OF NUCLEAR SYMMETRY ENERGY AROUND NORMAL DENSITY 

It is known that essentially all proposed nuclear interactions have been used in various many-body theories to 
predict the Esym{p) ISf. Instead of using pure model prediction, we use here the Egymipo) and 10 extracted from 
analyzing terrestrial nuclear laboratory experiments and astrophysical observations. Naturally, all analyses are based 
on some models and often different approaches are used in analyzing the same data or observations. For instance, at 
least 5 different models have been used to extract independently the Esymipo) and 10 from studying atomic masse. 
Remarkably, however, with very few exceptions, constraints on the Esym{po) and 10 from various analyses of the 



same or different experiments/observations overlaps closely. Listed in Table 1 are 24 sets of constraints including the 
4 astrophysical ones where only the upper or lower limit of 10 is given. We notice here that while some of the reported 
constraints provide both the upper and lower limits or the standard deviation together with the mean values, some 
do not provide any information about the associated uncertainties. This situation will be carried over into calculating 
the corresponding Usym{po,kp) and m^_p{pQ,S). We also notice that the constraints on the Esym{po) and 10 are 
treated here as quasi-data. Moreover, we use the most widely accepted empirical values for m^/jn and Ef{po) hi 
the formalism given in the previous section. Shown in Fig. 1 are the nucleon isovector potential (upper window) 
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FIG. 1: (Color online) Nucleon isovector potential Usymipojkp) (upper) and neutron-proton effective mass splitting 
mn-pipo,S)/5 (lower) at normal density of nuclear matter from 24 analyses of terrestrial nuclear laboratory experiments 
and astrophysical observations. 



and neutron-proton effective mass splitting (lower window) at po from the 24 constraints. We caution here that in 
cases where no error bar or range for the Esym{po) or 10 was given, only the error bar of the empirical value of 
rriQ/m is used in estimating the upper and lower limits of UsymipOikp) and ■m'^_ {pQ,6). In cases where only the 
upper or lower limits of 10 were given, the limiting values were indicated with arrows for comparisons. It is interesting 
to see that despite of the large uncertainty ranges of some of the constraints on Esymipo) and 10, the resulting 
mean values of Usyni{po,kp) and to*_ (po,<5) from different studies scatter very closely around their global averages 
of Usymipo, kp) =29 MeV and m*_p(po, S) = 0.24 • 6, respectively, indicating a high level of consistency of different 
studies. Moreover, the majority of the inferred m'^_p{po,6) are positive. While the mean values of Usymipo, kp) and 
iTT-n-pifOi ^) £ii'6 useful in their own rights, to use them reliably as a useful reference for calibrating nuclear many-body 
theories and much needed inputs for investigating many interesting issues in both nuclear physics and astrophysics. 



the community should strive at reducing the experimental error bars and providing more complete information about 
the associated uncertainties in extracting the Esym{po) and 10 using various models from the experimental data. In 
this regard, it is encouraging to note that some concerted efforts in this direction are under way by many people in 
both the experimental and theoretical communities investigating the density dependence of nuclear symmetry energy. 

IV. CONCLUSION 

Based on the Hugenholtz-Van Hove theorem, nuclear symmetry energy and the neutron-proton effective k-mass 
splitting are explicitly related to each other. Available constraints on the symmetry energy can be used to infer 
directly the poorly known but very important neutron-proton effective k-mass splitting in neutron-rich nucleonic 
matter. As an example, we have shown that the constraints on nuclear symmetry energy Esymipo) and its density 
slope 10 at pq from 24 recent studies of terrestrial nuclear laboratory experiments and astrophysical observations 
indicate consistently that the nuclear isovector potential and neutron-proton effective k-mass splitting at po are 
approximately Usymipo^kp) =29 McV and m'^_p{po,S) = 0.24 ■ 6, respectively. The need of reducing experimental 
error bars and proving more information about the associated uncertainties in extracting the Egymipo) and 10 are 
emphasized. 
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